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I. INTRODUCTION 



The process of photoassociation (PA) has served as a very useful tool to researchers 
in the field of ultracold Atomic, Molecular and Optical(AMO) physics.- The main advan- 
tage of this technique is its simplicity. The PA process provides an efficient, continuous, 
single-step method for the formation of alkali metal molecules in the ground X^S+ state 
and the metastable a'^S'^ state. ^ Detection of the molecules formed by PA and subsequent 
spontaneous emission is possible using several techniques, including resonance enhanced 
multiphoton ionization (REMPI). The REMPI process provides very rich spectra contain- 
ing information on not only the vibrational population distribution of the X state and the 
a state but also on the excited intermediate state via which ionization is performed. 

A major focus of the AMO community in recent years has been on the production of 
the absolute ground state of ultracold heteronuclear molecules. Such molecules may well be 
a stepping stone for the future realization of quantum computation and improved under- 
standing of many-body physics. Researchers all over the world have searched for different 
transfer pathways for the formation of these ground-state molecules, e.g. photoassociation^i^ 
and stimulated Raman adiabatic passage (STIRAP).-*^ 

Our group at the University of Connecticut has sought an efficient path for formation of 
absolute ground-state molecules of ^^K'^^Rb using single-step PA.™ Accurate knowledge of the 
excited electronic states is important to this research since the available ah initio potentials 
are not sufficiently accurate. Hence, we have performed high resolution spectroscopy of 
various electronic states (both excited and ground) of KRb over the years.->^"— In this 
paper, we report both experimental and theoretical investigations of the double-minimum 
3 electronic state of ^^K^^Rb, which may be useful in both the formation and detection 
of the X state and the a state. 



II. EXPERIMENT 

A detailed description of the experimental setup is given elsewhere.- Here we give a brief 
account. We start with a dual species "Dark-SPOT"— magneto-optical trap (MOT) of ^^K 
and ^^Rb. The average MOT density of ^^K is ~3xl0^° cm-^ and that of ^^Rb is -1x10^^ 
cm~^. After the two MOTs are carefully overlapped, we introduce a photoassociation laser 
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FIG. 1. Scheme for detection of the double-minimum 3 ^11 excited electronic state of KRb. The PA 
laser is detuned below the K(4s) + Rb(5pi/2) asymptote by 43.82 cm~^. The resulting a^S+ state 
molecules are then ionized via REMPI, where the 3^11 state is the resonant intermediate state. 

to form KRb molecules. The PA laser is a cw titanium: Sapphire ring laser pumped by a 
10 W 532 nm laser (Coherent Radiation model Verdi VIO). The typical output of the PA 
laser is ~1 W. The PA process produces excited state KRb molecules, which then undergo 
radiative decay to form molecules in the ground X^S+ state and the metastable a^E"*" 
state. Finally, the ground X and metastable a state KRb molecules are ionized using the 
REMPI method with time-of-fiight detection by a Channeltron charged particle detector. 
The ionization laser is a Continuum model ND6000 pulsed dye laser pumped at 532 nm by 
a doubled Nd:YAG laser. The REMPI laser has a linewidth of 0.5 cm~^ and produces 10 
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ns pulses of ~1 mJ energy, focused to a diameter (FWHM) of ~0.76 mm at the MOT. The 
KRb"*" ions thus produced are distinguished from other ions (Rb"*", K"*", Rb^) using time of 
flight mass spectroscopy. 

Electronic and vibrational spectroscopic information on the 3 ^11^ states is obtained by 
scanning the REMPI laser while keeping the PA laser fixed at the desired frequency. The 
experimental scheme is shown in Figure 1. The PA laser frequency is fixed at 12535.13 cm~^, 
populating the v' = 17, J' = 2 level of the 5(1) state, which at short range corresponds to 
the 2 state. For analysis of the 3 '^11^ potentials we use diabatic potentials obtained by 
modifying the existing ab initio potentials calculated by Rousseau et al}^ The lowest triplet- 
state potential correlating with two ground-state atoms is determined from experimental 
work.— 

III. MODIFIED AB INITIO POTENTIAL ENERGY CURVES 

Due to spin-orbit interactions, the 3 ^11 state splits into four separate potentials with 
= 0^,0", 1 and 2. We find that the f2 = O"*" and 0~ components are split by much less 
than the vibrational spacing, so it suffices to determine three potential curves for 3 ^Uq with 
i7 = 0, 1, and 2. Unfortunately, there are no direct ab initio calculations of the requisite 
Hund's case(a) potentials. Instead, the available ab initio potentials include a calculation of 
the potential at the short-internuclear range notation (appropriate for Hund's case (b)) for 
the 3^n state without spin-orbit terms—, as well as calculations of pure case (c) potentials 
that correspond to the = 0^ and 2 components, but not = 1. To construct approximate 
case (a) potentials for the analysis of our experimental data, we modify the existing Hund's 
case c) ab initio potentials where available, and use the case (b) potential to approximate 
the Hi component. 

Figure 2 shows the all of the available case (c) potentials converging to the K(4s)+Rb(4(i) 
and K(4s)+Rb(6s) atomic asymptotes^^, with panels (a), (b), and (d) depicting Q = 0+, 0^, 
and 2. There are numerous avoided crossings, and in the ) limit it is the diabatic 

continuation of each curve across these gaps that is needed. We have constructed approxi- 
mate potential curves for the Hq- , Ho+ , and H2 components of the 3 ^H state by using cubic 
splines to smoothly connect the adiabatic curves at the crossings. The dashed curves in Fig. 
2 show the resulting case (a) potentials. For the Hi component this procedure is impossible 
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FIG. 2. Plots of all available theoretical Hund's case (c) potentials arising from the K(4s)+Rb(4d) 
and K(4s)+Rb(6s) asymptotes . The = 0''^,0~,1 and 2 components are shown separately. In 
each case, the dashed line curve indicates the diabatic potential constructed to approximate the 
corresponding Vl component of the 3 '^11 state. The inset in (a) shows an expanded view of a typical 
avoided crossing between adiabatic curves. 

because the required 9(1) state has not been calculated, so instead we approximate the case 
(a) potential by directly using the calculated potential curve of the Hund's case (b) 3 '^11 
state, which is shown by the dashed curve in panel (c). This approximation is reasonable 
because the case (b) potential curve represents the weighted average of all spin-orbit com- 
ponents, which should lie close to the middle Hi component. It is also in good agreement 

5 




FIG. 3. Modified potential energy curves showing the spin-orbit components of the 3^11 state. 
The inset (a) shows the entire potential curves, while the main figure (b) shows an expanded view 
in which the $7 = 0, 1 and 2 components are distinct. The $7 = 0^ components, as expected, are 
nearly degenerate. 

with our experimental results, as we show in Section IV. 

In Figure 3, all four of the resulting approximate potential energy curves are shown in 
an expanded view, together with an inset showing their entire range. The near degeneracy 
of the f2 = 0"*^ components is evident, as is typical in Hund's case (a), and from now on 
we will refer to them together as i7 = 0. Our experimental observations confirm that the 
i7 = 0^ splitting is quite small. We have provided the listing of the modified potentials in 
the supplementary material.— 
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FIG. 4. Variations of the energy sphttings between different spin-orbit components of the 3^11 
state with respect to internuclear distance. The splittings are largest near ~5 A , corresponding 
to the bottom of the well. 



At long range the 3^11 potentials correlate to the K(4^Si/2)+Rb(4^Dj) atomic states, 
and the fine-structure splitting between J = 3/2 and J = 5/2 is only 0.44 cm~^. The much 
larger molecular fine-structure must gradually converge to this small value as the potentials 
approach their asymptotes. This is quantitatively represented in Figure 4, which shows that 
the spin-orbit splittings are largest (~40 cm^^) near the bottom of the well (~5 A) and 
decrease rapidly at long range. 
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IV. RESULTS AND ANALYSIS 



In this section, we first describe our PA and REMPI spectra, then analyze them using the 
potentials described in section III. We also discuss possible perturbations due to vibrational 
near-degeneracy between levels of the 3 ^IIj^ and the 4 states. 



A. Photoassociation Spectra 

Figure 5 shows PA spectra near 12535.1 cm~^. In these spectra two rotational bands 
are identified. The stronger KRb"*" signals correspond to J' = — 3 of the 3(0^) state with 
v' = 118 (the 1 '^n state at short range) and the weaker to J' = 2 — 4 of the 5(1) state with 
v' = 17 (the 2 state at short range), for which J' = 2 and 3 are most prominent. We were 
not able to observe the J' = 1 level of the 5(1) state in any of our PA scans. The set of three 
PA spectra with three different detection laser frequencies confirms the reproducibility of the 
spectral features both for the strong and the weak signals. The assignments of the PA spectra 
are based on previous experimental observation of the 3(0"'"), v' levels by photoassociation in 
a MOT.- The 5(1), v' levels were experimentally observed for higher values of J' by optical- 
optical double resonance polarization spectroscopy.-^ They were extrapolated to predict the 
positions of the 5(1), v' , J' = 1 — 4 levels.- 

The vibrational numbering of the 5(1) state is unambiguous since numerous levels were 
previously experimentally observed. However, for the 3(0"*") state, we have previously ob- 
served only five vibrational levels of the 3(0"*") state near the asymptote.- so the vibrational 
numbering is uncertain. Calculating vibrational levels of the 3(0"'") state using LEVEL^ 
and the ab initio potential, we find that v' = 118 of 3(0"*") is the closest match to the present 
photoassociation line. Thus we will refer to the 3(0"*") band reported in Figure 5 as f ' = 118, 
though the absolute numbering remains uncertain. 

In Appendix A we discuss possible problems due to overlaps with "hyperfine ghost" 
artifacts, arising due to a small population in the bright hyperfine states (F = 2 for '^^K and 
F = 3 for 85 Rb) of (4S) and ^^Rb (5S) in our dark SPOT, and we rule out the possibility. 
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FIG. 5. PA spectra indicating the overlapping 3(0"'") and 5(1) bands, with the REMPI detection 
laser frequency fixed at a) 16610.35 cm^^ b) 16624.76 cm^^ and c) 16644.01 cm~^. The three PA 
spectra confirm the reproducibility of the spectral features. The 5(1), J' = 2 line is from v' = 17 
and correlates to the short-range 2^11 state.— This is the excited state that produces a^S+ state 
molecules with v" = 15 — 25, and is marked by an asterisk. 



B. REMPI Spectra 

The REMPI spectra obtained by keeping the PA laser fixed at the 5(1), v = 17, J = 2 
level are shown in Figure 6. The spectra consist of transitions from the a^S+,f" = 15 — 25 
levels to the 3 ^Hq, = 2 - 12 levels, the 3 ^Hi, = 2 - 11 levels and the 3 ^Hs, i;' = 2 - 11 
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FIG. 6. REMPI spectra indicating spectral features corresponding to 3'^nQ-^a^S+ transitions 
(vertical lines) and atomic Rb nd^hs two-photon transitions. A variety of v" values from 14 to 
25 were observed, but for clarity only values from v" = 18 — 25 (above the noise level) are shown. 
The saturated atomic signal close to Rb 16d^5s corresponds to a ^^K atomic transition. The PA 
laser is fixed on the 5(1), v' = 17, J' = 2 level, at 12535.13 cm^^. 
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levels. The strongest transitions to the Q = 1 and 2 components of the 3 ^11 state originate 
from v" — 20 and 24 of the a state, and the strongest transitions to the Q = component 
of the S^n state originate from v" — 19 and 20. In Figure 6, we indicate only the transitions 
which are unambiguously assignable from the vibrational levels v" = 18 — 25 of the triplet 
metastable state to the spin-orbit components of the 3 ^11 state. The spectra are calibrated 
using atomic two-photon resonances, labeled by "no?", which are visible because the long- 
time tail of the atomic Rb+ time-of-flight distribution leaks into the time gate used for KRb+. 
The accuracy of the vibrational energies is determined by the REMPI laser hnewidth, which 
is 0.5 cm~^. The REMPI signal strength is dependent on the combination of the efficiency 
of a-state molecule formation in various v" levels by PA via the 5(1), v = 17, J = 2 level 
and the efficiency of the bound-bound transitions between the vibrational levels of the a 
and the 3 states. 



C. Analysis and Discussion 



TABLE I. Contributions of various J's (of S^Hq states) to the overall lineshape of the REMPI 
signals. 



J' 3^Uo 3^Ui 3^02 

10% 

1 20% 20% 

2 28.58% 40.48% 38.1% 

3 30% 21.67% 33.35% 

4 11.43% 17.8% 28.57% 



We have successfully assigned the vibrational levels = 2 — 11 of the 3^111 and 3^112 
state and f ' = 2 — 12 of the 3 "^Ilo state. From these observed transitions we obtain the 
total energy (E^iji) for a particular vibrational level v' which has a spread over its rotational 
levels (J') (these are not resolved due to the laser linewidth). The total energy is given by 
E„/ J/ = Tyi + J' [J' -\- 1)B„/ where T^ is the vibrational term energy and B„/ is the rotational 
constant. Hence to obtain T^ from E„/ j/ it is required to estimate the contributions of various 
J's to the overall lineshape of the transitions. In Table I, for each Q, component of the 3^11 
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FIG. 7. Comparison of theoretical and experimental values of the vibrational spacings, 
(cm^"'^) for the levels of S'^IIn states, where (a) shows the S'^Hq state, (b) the S^IIi state, and (c) 
the 3 ^112 state. The theoretical values are obtained using LEVEL^^ with the modified potentials 
reported in section III. The experimental values are obtained from the present work. The insets 
show an expanded view of the regions indicated by the rectangles (orange). 



state, we estimated contributions of various Js, using Honl-London factors.~ To obtain the 
numbers in Table I, we first calculated the spread of J's in the a-state resulting from the PA 
process and then the spread of J's in the 3 ^11 state due to the REMPI. From the numbers 
in Table I, it is possible to calculate the approximate shifts in E„/j/ values to obtain the T^ 
values. For the S'^IIo state the shift is ~8B^/, for the S^IIi state the shift is ~8.99B„/ and 
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for the 3^112 state the shift is ~12Bt,/. For example, the theoretical B^/ values (calculated 
using LEVEU^) for v' = 3 of 3^00, 3^Ui and 3=^n2 states are 0.02235 cm-\ 0.02309 cm'^ 
and 0.02652 cm~^ respectively. Then the required shift in Ei,/=3 j/ to obtain Tv>=3 for the 
3^no, 3 ■^IIi and 3^112 states are 0.18 cm^^, 0.21 cm~^ and 0.32 cm~^ respectively. For our 
experiment this shift is well within the laser linewdith, so there is no significant difference 
between Et,/j/ and T^. 

TABLE II. Comparison of experimental and theoretical term energies T^' values of the vibrational 
levels = — 12 of the S'^IIq- state. We also compare the theoretical vibrational term energies of 
the S'^IIg- and 3 '^110+ state, for the vibrational levels v' = — 12. The term energies are reported 
with respect to the K(4Si/2)+R'b(5Si/2) asymptote and are in units of cm^"*^. 



v' 


Experiment 


Theory 














rr-i 















16600.75 


38.32 


0.49 




1 






16639.07 


36.77 


-0.56 




2 


16541.96 


33.75 


16675.83 


33.89 


-1.56 


133.88 


3 


16575.70 


25.73 


16709.73 


24.85 


-0.07 


134.02 


4 


16601.44 


12.62 


16734.58 


13.48 


-0.60 


133.15 


5 


16614.06 


23.74 


16748.07 


23.59 


-0.46 


134.01 


6 


16637.79 


23.49 


16771.66 


23.38 


-0.62 


133.86 


7 


16661.28 


25.04 


16795.04 


24.98 


-1.02 


133.75 
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16686.33 


27.78 


16820.01 


26.22 


-0.85 


133.68 
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16714.11 


25.40 


16846.23 


26.99 


-0.69 


132.13 


10 


16739.51 


28.23 


16873.22 


28.18 


-0.28 


133.71 


11 


16767.74 


29.31 


16901.40 


29.36 


0.10 


133.66 


12 


16797.05 




16930.76 


30.04 


0.66 


133.71 



^ AT„/ = T„,{0-) - T^,{0+) [both theoretical values; 
^ AT„/ = T^,/ (Theory) -T^,.(Expt.) 



The vibrational term energies (T„/) and the vibrational spacings (AG„+i/2) of each of 
the Q components of the 3 ^11 states are tabulated in Tables II, III and IV. In these tables 
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TABLE III. Comparison of experimental and theoretical term energies T^i of the vibrational 
levels v' = — 11 of the S^IIi state. The term energies are reported with respect to the 
K(4Si/2)+Rb(5Si/2) asymptote and are in units of cm~^. 
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10 
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we also compare the experimental values of T^' and AGi,+i/2 with theoretical predictions 
(calculated using LEVEL^) based on the approximate curves described in section III. The 
calculated values of Ty/, are higher than the experimental values obtained from experimental 
data by an average of ~133.6 cm^^ for the S^IIo state, ~139.3 cm~^ for the S^IIi state and 
~138.1 cm^^ for the 3 '^112 state. However the experimental vibrational spacings, AG^,+i/2, 
are within 6 — 8% of the theoretical values, as can be seen in Figure 7. The dip in the plot 
of AG^+1/2 vs. energy near 16600 cm~^ (Figure 7) is due to the double-minimum character 
of the potential, evident in Figure 3. 

We conclude that the theoretical potentials reported in section III when shifted down by 
134-139 cm~^, provide an accurate potential energy curve for the double minimum 3^11^ 
state in the region of the observed vibrational levels. In Figure 8 we show the 3 ^IIi state, 
where the dashed curve indicates the ab initio potential and the solid curve shows this 
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TABLE IV. Comparison of experimental and ab initio term energies Tyi of the vibrational levels 
= — 11 of the 3 '^112 double minimum state. The term energies are reported with respect to 
the K(4Si/2)+Rb(5Si/2) asymptote and are in units of cm"-*^. 
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potential shifted down by 139 cm . 

Also shown in Figure 8 is the portion of the 4 ^S"*" potential curve that overlaps with the 
3^nf7=i state. The 4^S+,t/ = 1 — 16 levels were previously experimentally observed^°. As 
can be seen in the figure, some vibrational levels of these two states lie in close proximity. 
Pairs closer than ~2 cm~^ include the 4^S"'",f' = 15 — 3^11, f' = 6 and the 4^S+,f' = 
16 — 3 '^n, v' = 8. This suggests that these levels may appreciably perturb each other. 
Perturbations may also be significant for the unobserved f' = 17 level of the 4^S+ state 
and the v' = 10 level of the 3 ^11 state. However, these perturbations are not obvious in the 
data in Table III given the ~0.5 cm~^ uncertainty from the pulsed laser linewidth. Higher 
resolution studies would be desirable. 

We have not observed levels above v' = 12 of the 3 ^Hq states, only because REMPI scans 
were not performed beyond the reported range. Our inability to observe the levels v' = 
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R(A) 

FIG. 8. The solid (black) double minimum potential curve shows the S^IIi state, shifted down 
by 139 cm^^ (indicated by the double headed arrow). The original approximate potential curve is 
indicated by the dashed (red) line. The solid horizontal lines indicate the vibrational level positions 
of the S^Ili state for which v' = and 1 levels are calculated and not experimentally observed. 
The other solid (blue) potential curve is the 4^S"*' state. The horizontal dashed lines indicate the 
vibrational levels of the 4^S^ state, that have been experimentally observed^. Note the near- 
degeneracy of the v' = 6 and 8 levels of the S^lli state with the v' = 15 and 16 levels of the 4^S"'" 
state, respectively. 

and 1 is probably due to the lack of Franck- Condon overlap of these two levels with the 
a^Il'^,v" levels. This is evident if we focus on the inner turning points (R^-) of the a^S"*" 
and 3 ^11 states. As previously discussed, the a ^E"*" state is mainly populated in the range 
v" = 15 — 25. The corresponding inner turning points are Ri5//_=4.99 A and R25"_=4.92 
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FIG. 9. The FCFs for S^Hi,?;' = - ll^a'^T.+,v" = 20 and 24, calculated using LEVEL^^. 

A. However Figure 8 indicates that the S'^HijV' = and 1 levels lie deep in the outer well 
at a distance larger than 4.99 A, yielding a poor Franck-Condon overlap. These Franck- 
Condon Factors (FCFs) have been calculated using LEVEL^ and are shown in Figure 9 for 
transitions from a ^5]+, v" = 20 and 24 to 3 ^Hi, t;' = - 11. The FCFs for the 3^Ui,v' = 
and 1 levels are negligible and for the v' = 2 level they are relatively weak, in agreement 
with our experimental observations and the above argument. 

Also, to help characterize this double-well potential. Figure 10 shows the vibrational 
wavefunctions of this state for all the observed vibrational levels. The outer well is deeper 
than the inner well as shown in Figure 8. Thus the f ' = — 2 levels belong entirely to the 
outer well, while the v' = 3 wavefunction has a small penetration into the inner well, which 
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FIG. 10. The vibrational wavefunctions of the S^Ili double-minimum state for v' = — 11, 
calculated using LEVELS^. 

is apparent from both its level position (Figure 7) and wavefunction (Figure 10). Above 
v' = 3, the amplitude of the wavefunction gradually increases in the inner well region and 
from v' = 7 onwards the wavefunctions have a smooth envelope spanning both wells. 

The other spin-orbit components (f2 = and 2) of the 3 ^11 state should follow very 
similar arguments to those just discussed for f2 = 1. 
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V. CONCLUSION 



In conclusion, we have used REMPI of ultracold '^^K^^Rb to observe the double minimum 
3 ^Ilf^ states vibrational levels f ' = 2 — 1 1 of the f2 = 1 and 2 components and f ' = 2 — 12 
of the = spin component. Molecules in the a^S+ state are formed by use of a PA 
laser fixed at 12535.13 cm~^ which corresponds to the 5(1), f' = 17, J' = 2 level. We are 
able to assign the spin-orbit components of the 3 ^11 state based on theoretical potentials 
obtained by modifying the existing ab initio potentials.— Also, we predict that the 
and the 3 ^Uq states may perturb each other at certain vibrational levels which needs further 
investigation. 

This new spectroscopic information for the excited electronic states of KRb provides 
better understanding of this molecule. It removes the ambiguity in the electronic and vibra- 
tional energies of the excited states in this energy region, thus providing new opportunities 
for different transfer pathways for future experiments on ultracold molecules. 
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Appendix A 

The 5(1), v' = 17 rotational band lies very close to the predicted positions of '^^K hyperfine 
ghosts of 3(0'''), f' = 118. Here, we provide further analysis to confirm our assignments of 
the 5(1), = 17, J' = 2 - 4 levels. 

Figure 11 is the PA spectra shown in Figure 5b with additional labels showing the pre- 
dicted positions for '^^K hyperfine ghosts of 3{0'^),v' = 118. As can be seen, the 5(1), J' = 2 
level lies just below the hyperfine ghost position for 3{0^),J' = 1, while the 5(1), J' = 3 
is just above the hyperfine ghost position for 3(0"'"), J' = 3. Furthermore, there is no clear 
evidence of any hyperfine ghosts for 3(0"*"), J' = 0, 2, where no near overlaps occur. 

To investigate further, we acquired a pair of REMPI spectra over the same wavelength 
region, with PA first to 5(1), J' = 2 and then to 3(0^), J' = 1. Figure 12 shows the resulting 
spectra. If one of the two lines is the hyperfine ghost of the other, then their REMPI spectra 
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FIG. 11. Redrawing of the PA spectra shown in Figure 5b with additional dotted lines indicating 
the predicted positions of ^^K hyperfine ghosts of the 3(0^), v' = 118 rotational band. 



should be almost identical. However, Figure 12 clearly shows that the two REMPI spectra 
are quite different. 

Hence, we conclude that our assignments of the rotational band corresponding to 5(1), 
v' = 17 (refer to Figure 5) is correct and they are not hyperfine ghosts of the rotational 
bands of 3(0+), = 118. 
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16580 16600 16620 16640 

REMPI Laser Wavenumber (cm"^) 

FIG. 12. Comparison of the REMPI spectra with PA fixed at 5(1), J' = 2 (black, narrower width) 
and at 3(0+), J' = 1 (red, thicker width). The "nd" notations indicates the atomic two-photon 
transitions that leaked into the molecular channel. These lines serve as precise internal references 
for calibration of each spectrum. 
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